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in: 

We have calculated the ratios R {P ) = T(P -> eV e \y\)/T(P -> £lV M [y]) (P = 7t,K) in Chiral Per- 
turbation Theory up to G{e 2 p% finding = (1.2352 ±0.0001) x 10~ 4 and = (2.477 ± 
^T* ■ 0.001) x 10~ 5 . This observable is helicity suppressed in the Standard Model, so that it is a sen- 

£f~) . sitive probe of all Standard Model extensions that induce pseudoscalar currents and nonuniversal 
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corrections to the lepton couplings. Ongoing experimental searches plan to reach uncertainties 



that are comparable to these results. At the moment is in agreement with the final result 
by the KLOE Collaboration at DAFNE and it is at 1 .4 <T of the preliminary result by the NA62 

(k) 

Experiment at CERN. New measurements of /T^ are under way by the PEN Collaboration at 
PSI and by the PIENU Collaboration at TRIUMF. 
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1. Motivation 

The ratio R^) = T(P -> ev e [y])/T(P jUV M [y]) (P = 7t,K) of leptonic decay rates of light 
pseudoscalar mesons is helicity-suppressed in the Standard Model, due to the V —A charged current 
coupling. It is therefore a sensitive probe of all Standard Model extensions that induce pseudoscalar 
currents and non-universal corrections to the lepton couplings [j|]. Attention to these process has 
been payed in the context of the Minimal Supersymmetric Standard Model, with [Q] and without [Q] 
lepton-flavor-violating effects. In general, effects from weak-scale new physics are expected in 
the range (AR e /n)/R e /n ~ 10~ 4 — 10~ 2 and there is a realistic chance to detect or constrain them 
because of the following circumstances: 

(k) (71) 

i) Ongoing experimental searches plan to reach a fractional uncertainty of (A/?^)//?^ ~ 
5 x 1(T 4 [g |] and (AR^h/R^l ~ 3 x 1CT 3 [@, 0], which represent respectively a factor of 



e/flJ/ e/fl 

around 5 and 10 improvement over former errors [Jsj, |9[ [R], O, 12, 13] 



ii) At the same time, the Standard Model theoretical uncertainty can be pushed below this 

level, since to a first approximation the strong interaction dynamics cancels out in the ra- 
ff) 

tio 7?V and hadronic structure dependence appears only through electroweak corrections. 
Indeed, the most recent theoretical predictions read R^) = (1.2352 ±0.0005) x 10~ 4 [0], 



R 1 ?) = (1.2354 ±0.0002) x 10~ 4 @, and R { * ] = (2.472±0.001) x 10~ 5 InRef. [0] 



a general parameterization of the hadronic effects is given, with an estimate of the leading 
model-independent contributions based on current algebra [|l6|]. The dominant hadronic un- 
certainty is roughly estimated via dimensional analysis. In Ref. [O], on the other hand, the 
hadronic component is calculated by modeling the low- and intermediate-momentum region 
of the loops involving virtual photons. 

2. The Standard Model prediction 



1— -1 I — 1 (p\ 1 — 1 

In Refs. |[17|, |18|] we have analyzed /?• / within Chiral Perturbation Theory []19D, the low- 



energy effective field theory of QCD. The key feature of this framework is that it provides a con- 
trolled expansion of the amplitudes in terms of the masses of pseudoscalar mesons and charged 
leptons (p ~ m K ^ K t/A x , with A^ ~ 4nF n ~ 1.2GeV), and the electromagnetic coupling (e). Elec- 



tromagnetic corrections to (semi)-leptonic decays of K and n have been worked out to 6[e 2 p 2 ) [ ]20| , 

p) 

/n 



21], but had never been pushed to &(e 2 p A ), as required for R^f) in order to match the experimental 



accuracy. 

Within the chiral power counting, R e ^ is written as: 



Jf(P) _ j(0),(f) /, , A \ 



l+A^+A^ + A^U-. 

e z p z e L p* e^p a 



(2.1) 



being R^ the well known tree-level expression: 
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(P=n) 


(P = K) 







(7.84±0.07 r ) x 10~ 2 


c (P) 


5.2±0.4 L ,±0.01 r 


4.3±0.4 L ,±0.01 r 


c (P) 


-10.5 ±2.3 m ±0.53^ 


-4.73±2.3 ra ±0.28 L9 


4 ho 


1.69 ±0.07^ 


0.22±0.01 L9 



(p) (p) 

Table 1: Numerical values for c 2 3 4 and c 2 , for P — 7t,K. The uncertainties correspond to the input values 
L r g (p = m p ) = (6.9 ±0.7) x 1CT 3 , 7= 0.465 ±0.005 [Q, and to the matching procedure (m), affecting only 

JP) 



At the level of uncertainty considered here, one needs to include higher order long distance correc- 



tions | ]14p and their effect amounts to the factor 1 +An in ( |2. 1[ ), 



1+4 



m e I 



LL 



l_^log^ii 



1.00055 



(2.3) 



The leading electromagnetic correction A^ 2 corresponds to the point-like approximation for pion 
and kaon, and its expression is also well known jl4] , |20| , ^2| ] : 



A (P) 

e z p- 



mz 



mi 



13 



m. 
> P 

"mp 
\9z 



8-5z 



8(1 -z) 4(1 - z ) 



zlogz 



2 i±£L/ 2 (l- z ) . 



2± -^logz ) log(l -z) 



(2.4) 



(2.5) 



(p) 

The structure dependent effects are all contained in A^ 2 ; 4 and higher order terms, which are the main 



subject of Refs. [jn, 18]. Neglecting terms of order (m e /m p ) 2 , the most general parameterization 
of the next-to-leading chiral contribution can be written in the form 



,(P) 



^ ^log^+cr+cfK/^) +-^4 p 

K ml \ z ml 1 / 71 ml z 



log 



m^ 

ml 



(2-6) 



(p) 

which highlights the dependence on lepton masses. The dimensionless constants c\ 3 do not de 

(p) 

pend on the lepton mass but depend logarithmically on hadronic masses, while c\ (m^/mp) 







in) 

as win — » 0. (Note that our do not coincide with C23 of Ref. Q14|], because their C3 is not 



constrained to be mi -independent and contains in general logarithms of rri£.) 

Let us note that the results for c^l^ and cip depend on the definition of the inclusive rate 
r(P — > -fv^[y]). The radiative amplitude is the sum of the inner bremsstrahlung (Tib) compo- 
nent of @{ep) and a structure dependent (Tsd) component of &(ep 3 ) p3[]. The experimental 



(71) 1 1 1 

definition of R / is fully inclusive on the radiative mode, so that A V 2 ; 4 receives a contribution 



e/n 



from the interference of Tjb and Tsd- Moreover, in this case one also has to include the ef- 



fect of A 



(*) 



e 2 p6 x \Tsd\ , that is formally of 



e p ), but is not helicity suppressed and behaves 



3 
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(P=n) 


(P = K) 


41 (*) 


-3.929 


-3.786 


41 (*) 


0.053 ±0.011 


0.135±0.011 


Ag, (%) 


0.073 










All (%) 


0.055 


0.055 



Table 2: Numerical summary of various electroweak corrections to R),^- The uncertainty in A^ 2 4 corre- 
sponds to the matching procedure. 



as A e 2 p 6 ~ a / K {mp / My) A (mp/m e ) 2 . On the other hand, the usual experimental definition of 

(K) 

is not fully inclusive on the radiative mode. It corresponds to including the effect of Tjb in A , 2 

€ p 

(dominated by soft photons) and excluding altogether the effect of Tsd- consequently c„ ^ c„ . 

(p) (p) 

The expressions of the constants c 2 3 4 and cj, are shown in Refs. §\J\, |T^] and their numerical 

values are reported in table []]. Note that to this order in Chiral Perturbation Theory, Rgil features 
both model independent double chiral logarithms (previously neglected) and an a priori unknown 

low-energy constant. By including the finite loop effects and estimating the low-energy constant 

(p) 

via a matching calculation in large-Afc QCD, we thus provide the first complete result of /?V to 
G(e 2 p A ~) in the effective power counting. Most importantly, the matching calculation allows us to 
further reduce the theoretical uncertainty and put it on more solid ground. 

In table | we summarize the various electroweak corrections to R^J'^ ■ Applying these we 
arrive to our final results: 

R^l = (1.2352 ±0.0001) x 10~ 4 , (2.7) 



e/n 

R<& = (2.477 ±0.001) x 10~ 5 . (2.8) 



<7l) 

The uncertainty we quote for R V is entirely induced by our matching procedure. However, in the 

( K) 

case of R e /L w e have inflated the nominal uncertainty arising from matching by a factor of four, to 
account for higher order chiral corrections, that are expected to scale as A^ 4 x m\j (4nF) 2 . 



3. Discussion 



3.1 Comparison to previous theoretical predictions 



Our results have to be compared with the previous theoretical predicitions of Refs. []14j] and 
[15], which we report in table ||. 

i) R\ J is in good agreement with both previous results. 

ii) There is a discrepancy in Zr ,^ that goes well outside the estimated theoretical uncertainties. 
We have traced back this difference to two problematic aspects of Ref. [[1^]. The leading 
log correction All is included with the wrong sign: this accounts for half of the discrepancy. 
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10*.*$ 


10= -<fj 


This v 
Ref. I 
Ref. [ 


/or 
14 
15 


k 

1 

] 


1.2352 ±0.0001 
1.2352 ±0.0005 
1.2354 ±0.0002 


2.477 ±0.001 
2.472 ±0.001 



Table 3: Comparison of our result with previous theoretical predictions of R 



(P) 
•III- 



The remaining effect is due to the difference in the next-to-leading order virtual correction, 



for which Finkemeier finds A 



(K) 



0.058%. We have serious doubts on the reliability of 



this number because the hadronic form factors modeled in Ref. [15] do not satisfy the cor 



rect QCD short-distance behavior. At high momentum they fall off faster than the QCD 

(K) 

requirement, thus leading to a smaller value of A 2 4 compared to our work. 

e p 



3.2 Comparison to experiments 



3.2.1 * 



(*) 

e/n 



The three most recent measurements of R , are mutually consistent 



e/n 

*$|Bryman = (1.218 ± 0.014) X 10~ 4 Ref. § , 

*i/ilBritton = (1.2265 ±0.0034 stat ±0.0044 syst ) x 10~ 4 
*^|czapek = (1.2346 ±0.0035 stat ±0.0036 syst ) x 10~ 4 



(1.227 ±0.006) x 10~ 4 
(1.235 ±0.005) x 10~ 4 



Ref. @ , 

Ref. . 

(3.1) 



These measurements are in agreement with the theoretical prediction and they give the PDG av- 
erage (1.230±0.004) x 10~ 4 pBj]. Note that it is less accurate than the prediction by a factor of 
around 40. As it has been indicated previously, the experiments ruled by the PEN Collaboration at 
the Paul Scherrer Institute [Q] and by the PIENU Collaboration at TRIUMF [||] are under way and 
are expected to improve significantly the uncertainty. 



3.2.2 * 



(K) 



The old measurements of *^P in the seventies are also mutually consistent and in agreement 
with the Standard Model prediction, 



*i/ilciark = (2.42±0.42)xKT 5 

*S|Heard = (2.37 ± 0. 17) X 10 
*^| H emtze = (2.51 ± 0.15) X 10 



-5 



-5 



Ref. , 
Ref. [|E 
Ref. C- 



(3.2) 



Using these data the world average reads (2.45 ± 0. 1 1) x 10 [25]. Again it is much less accurate 
than our theoretical prediction, by a factor of around 100. Recent experiments have been improving 
significantly the uncertainty: 
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Figure 1: Summary of r V measurements (picture from 



i) The KLOE Collaboration has recently published the measurement performed at DAFNE [0], 
R§1 Ikloe = (2.493 ± 0.025 stat ± 0.019 syst ) x lO" 5 = (2.49 ± 0.03) x 10" 5 , (3.3) 



in agreement with our Standard Model prediction of 2.8. The achieved precision is 1.2%, 



improving the precision of the former world average by a factor of 4. 

ii) On the other hand the NA62 experiment at CERN has also recently announced its preliminary 
result [0], 

#2|na62 = (2.500 ±0.012 sta t± 0.01 l syst ) x 10~ 5 = (2.500 ±0.016) x 10~ 5 . (3.4) 

The uncertainty is now 0.64%. Note that the whole 2007-08 data sample is supposed to allow 
pushing the uncertainty down to 0.4%. This result is compatible with the KLOE one and it 
is also in agreement with the theoretical prediction of ^ (at 1.4 a). 

With these new results, and until the final result of NA62 arrives, using a simple weighted mean 
the new world average reads 

|wa = (2.498 ± 0.014) x 10" 5 , (3.5) 

which is in agreement with the Standard Model result (at 1 .5 a). A summary of R^ measurements 
is presented in figure [ij 
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